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Abstract 

^ , We study the distribution of color electric flux of the three-quark system in quenched and full 

^ ' QCD (with Nf = 2 flavors of dynamical quarks) at zero and finite temperature. To reduce ultra- 

violet fluctuations, the calculations are done in the abelian projected theory fixed to the maximally 
abelian gauge. In the confined phase we find clear evidence for a Y-shape flux tube surrounded and 
formed by the solenoidal monopole current, in accordance with the dual superconductor picture 
of confinement. In the deconfined, high temperature phase monopoles cease to condense, and the 
distribution of the color electric field becomes Coulomb-like. 

PACS numbers: 11.15.Ha, 12.38.Aw, 12.38. Gc 
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I. INTRODUCTION 



So far most investigations of the static potential, and the dynamics that drives it, have 
concentrated on the quark-antiquark (QQ) system, while little is known about the forces 
of the three-quark (3Q) ensemble. For understanding the structure of baryons and, in 
particular, for modelling the nucleon it is important to learn about the forces and the 
distribution of color electric flux in the 3Q system as well. A particularly interesting question 
is whether a genuine three-body force exists and the confining flux tube is of Y-shape, or 
whether the long-range potential is simply the sum of two-body potentials, in agreement 
with a A-Ansatz, resulting in a flux tube of A-shape. By a flux tube of Y- and A-shape 
we understand a flux tube between the three quarks having shortest possible length and a 
junction, and a flux tube constructed out of three quark-antiquark flux tubes taken with a 
factor |. 

Several lattice quenched QCD studies report evidence for a A-type long-range poten- 
tial , while others claim a genuine three-body force 0,0]. In Ref. jj^ various patterns of 
the three-quark system were considered with the distance between quarks in an equilateral 
triangle, d, up to 0.8 fm. It was found that at large distances the y-Ansatz gives a better 
description of the three-quark potential than the A-Ansatz. On the other hand, the authors 
of Ref. Pi found that at distances d < 0.7 fm the three-quark potential is described quite 
well by A-Ansatz, while it rises like the y-Ansatz at larger distances, 0.7 < d < 1.5 fm. 

The F-Ansatz is also being supported by the field correlator method The difference 
between a A- and Y-shape potential is rather small and difficult to detect, because the 
underlying Wilson loop decays approximately exponentially with increasing interquark dis- 
tance. A computation of the distribution of the color electric flux inside the baryon might 
help to resolve this problem. 

In this paper we shall study the static potential and the flux tube of the system. 
The long-distance physics appears to be predominantly abelian - being the result of a yet 
unresolved mechanism - and driven by monopole condensation. The use of abelian variables 
is an essential ingredient in our work, as it leads to a substantial reduction of the statistical 
„oise. Preliminary results of this ..vestigat.ou Kave bee. reported „r Ref. Q. 

The paper is organized as follows. In Section 2 we describe the details of our simula- 
tion, including the correlation functions that we are going to compute. The results of the 
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calculation are presented in Sections 3 and 4. Section 3 is devoted to the study of the 3Q 
system at zero temperature, while Section 4 deals with the finite temperature case. Finally, 
in Section 5 we conclude. 



II. SIMULATION DETAILS 

We employ the Wilson gauge field action throughout this paper. In our studies of full 
QCD we are using non-perturbatively 0(a) improved Wilson fermions, 

Sf = S'p^ - ^ng cswa^ ^ ^{s)crf,^F^^{s)ilj{s) , (1) 

s 

with Nf = 2 fiavors of dynamical quarks, where Sp'^ is the ordinary Wilson fermion action. 
Further details of the dynamical runs are given in 

The system of three static quarks propagating from A to B may be described by the 
'baryonic' Wilson loop 

= —^£ijk£i'jik'Uiii{Ci)Ujji{C2)Ukk'iC-i) 1 (2) 

where 

U{C)=\{U{s,^) (3) 

is the ordered product of link matrices U G SU (3) along the path C, as shown in Figjn] The 
potential energy of this system is given by 

V = -^ Mm log(W3Q), (4) 

Lj- being the temporal extent of the loop. 

In the following we shall concentrate on abelian variables, referring to the maximally 
abelian gauge (MAG), and being obtained by standard abelian projection 12 1. To fix 
the MAG jlj], we use a simulated annealing algorithm described in [l^. We write the 
abelian link variables as 

u{s,fi) = diag(Mi(s,/i),M2(s,/i),M3(s,/U)) , 

(5) 
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FIG. 1: Three quark Wilson loop. 



with 

1 ^ 

9i{s,n) = arg(t/ii(s,/i)) - - ^ arg(%(s, /x)) 



mod 2iT ' 



(6) 

4 4 

^i(s,/i) G [--TT, -Tt] . 

They take values in U{1) x and under a general gauge transformation they transform 

as 

m(s, fx) — i> (i(s)^'u(s, /i)(i(s + jl) , 
(i(s) = diag(exp(i ai(s)),exp(i 0:2(5)), (7) 
exp(-i(ai(s) + 02(5)))) . 

The abelian Wilson loop is given by 

WsQ = ^ h,k\u,{C,)u,{C2)uk{Cs) , (8) 

where u{C) is the abelian counterpart of W^q is invariant under gauge transformations 
0. 

The physical properties of the 3Q system can be infered from correlation functions of 
appropriate operators with the corresponding Wilson loop. Abelian operators take the form 



0{s) = diag(Oi(s), 02{s), 0^{s)) e U{1) X U{1) . (9) 
For C-parity even operators O, like the action and monopole densities, the correlator of 



0{s) with the abelian Wilson loop is given by tI. I14 | 



iOis)hQ = -^y^F^ - (^) • (10) 

For C-parity odd operators, like the electric field and monopole current which carry a color 
index, the correlator is defined by 

{0{s)hQ-^_ , (11) 

where summation over i,j,k is assumed. It is natural to use Wilson loop to study the 
static potential at zero temperature since it gives directly a singlet potential. The Polyakov 
loop correlator gives in general a color-averaged potential, i.e. a mixture of the singlet and 
octet potentials, see e.g. P|. At nonzero temperature one can use only the Polyakov loop 
correlator to study the static potential and we use the product P^q of three Polyakov loops 
closed around the boundary as baryonic source instead of W^q: 

PSQ = ^ \^^Jk\ (^^iSl) £,(82) 4(^3) , (12) 

where 

^.(^1= Y[u,{s,s,,4) (13) 

S4 = l 

is the abelian Polyakov loop, Lt being the temporal extent of the lattice here. The correlators 
of 0{s) with P3Q are defined analogously to (fTUj) and (fTT|) . 

The observables we shall study are the action density , the color electric field E^^ and 
the monopole current k^'^. The action density is given by 

PA^(^) = f E<^°<^^(^'/^''^)))3Q' (14) 

where 

Oi{s,p,u) = aTg{ui{s,n,u)) , ^^^^ 
Ui{s, fi, v) = Ui{s, fi)ui{s + fi, i^)ul{s + 0, fi)ul{s, v) 
is the plaquette angle. The color electric field and monopole current correlators are defined 
by 

E3«(s, i) = i (diag(^i(s, 4, z), ^2(5, 4, z), ^3(s, 4, i)));q (16) 



and 

/i) = 27ri (diag(A;i(*s, fi), k2{*s, fesf fi)))sQ , (17) 
where ki{*s,fi) is the monopole current jiol . 14 1. 



The calculations in full QCD at zero temperature are performed on the 24^ 48 lattice at 
(3 = 5.29, K, = 0.1355, which corresponds to a pion mass of m.^/mp ^ 0.7 and a lattice 
spacing of a/ro = 0.18 (i.e. a = 0.09 fm assuming tq = 0.5 fm). The calculations in full 
QCD at finite temperature T are done on the 16'^ 8 lattice at /? = 5.2 for various hop ping 
parameters ranging from k = 0.1330 to k = 0.1360, which covers the temperature range [l^ 
0.8 < T/Tc < 1.2. The critical temperature Tc corresponds to k = 0.1344(1). At this k we 
find m-j^/mp ~ 0.77. For comparison, we also did quenched simulations at zero temperature 
on the 16'^ 32 lattice at /3 = 6.0. At this (3 the lattice spacing is a/ro = 0.186, i.e. it is 
roughly the same as on our full QCD lattices. To reduce the statistical noise we smeared the 
spatial links of the abelian Wilson loop using 10 sweeps of APE smearing with a = 2, 
where a is a coefficient multiplying the sum of staples. 



III. STATIC POTENTIAL AND BARYONIC FLUX AT ZERO TEMPERATURE 



The minimal Y-type distance between the quarks, i.e. the sum of the distance from the 
quarks to the Fermat point is ^ 



1/2 



(18) 



where marks the position of the i quark, Vij = \ri — rj \ and Sa is the area of the triangle 
spanned by the three quarks. The y-Ansatz predicts that the confining part of the baryonic 
potential is a'^Ly, with string tension a'^ equal to the QQ string tension 17 1: 



a 



3Q _ ^QQ 



Y 



a 



The full expression describing both large and small distances is 



Y 



T/3Q 



i<j 



3Q 



(19) 



(20) 



where, similarly to the QQ static potential, V^^ is a selfenergy term, the Coulomb term 
with effective coupling a^^ comprises one gluon exchange as well as a Liischer term, recently 



derived for the baryonic string in [18|], and the confining term has string tension a'^ . The 
A-Ansatz prediction jl^ is that the confining part of the potential is proportional to the 
perimeter of the triangle formed by the quarks 

L^ = J2\^^-f,\. (21) 

with string tension 

af = la^^. (22) 

The short distance part is of the same form as in eg. ()2()|1 . Thus the full expression for the 
A-Ansatz potential is 

V'^^iL^) = V^^-Y, — + <yfL^ ■ (23) 

For short distances perturbation theory arguments relate the selfenergy and the Coulomb 
term coefficient to those of the QQ static potential j^: 

y^Q^h^QQ^ ^3Q^1^QQ_ (24) 

On the other hand fitting the numerical data including both long and short distances by 
(I2()|) or by ()23|1 one may find results which differ from ()24|1 . e.g. due to the Liischer term 
contribution. In Ref.^] a rough agreement between the fit parameters and ()24|) has been 
found for both F-Ansatz and A-Ansatz fits. 

In Fig 121 we show the baryon potential as a function of Ly. An unphysical constant 
has been subtracted from the potentials. For equal distances between the quarks, i.e. 



r,; - TA 



d = for Vi 7^ j, eq. piHl becomes 



y3Q(LY) = K,"^ - SVS^ + Ly . (25) 

Ly 

Fitting our data for three quarks in equilateral triangles by eq. ()25|) for distances d < 0.75 fm 
we found the abelian string tension 0"^*^^^ = 0.038(1) and 0.0395(12) for the quenched and 
full theory, respectively. These values agree within error bars with the abelian string tension 
for the QQ fiux tube a^^ = 0.039(1) and 0.0402(11) Q thus supporting the F-Ansatz. 
Note, that both 3Q and QQ abelian string tensions are slightly higher in full QCD. We found 
values for the selfenergy and the Coulomb term coefficient smaller than prescribed by ()24j] : 
Y^Q/yQQ = 1.28(5), = o.27(4) in full QCD and V^'^ /Vo^"^ = 1.31(6), a^Q/^QQ = 

0.31(4) in quenched QCD. The fits are also shown in Fig|21 



Full H 




1 2 



FIG. 2: The abelian baryon potential in full and quenched QCD. 
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FIG. 3: Comparison of the abelian and SU(3) baryon potential in the quenched approximation on 
the 16^ 32 lattice at /? = 6.0. The SU(3) potential is taken from J|. 



In FigO the abelian and the nonabelian quenched potentials are plotted together with 
respective fits. The data for the nonabelian potential is taken from Ref. ^. Comparison of 
"^yab with the SU(3) result ^ gives <yY^g,\^/<y'^ = 0.83(3), which lends further support to the 
hypothesis of abelian dominance. 

If the confining flux is of Y-shape we would expect the long-distance part of the potential 
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FIG. 4: The abelian baryon potential in full QCD, together with its monopole and photon part, 
as a function of Ly (top) and La (bottom), respectively. The curves show fits (top) and 
(bottom) to the data with equal distances between the static sources. 

to be a universal function of Ly. In Fig0]we plot the abelian potential as a function of Ly 
(top), and as a function of La (bottom). The data show a universal behavior when plotted 
against Ly. This is to a lesser extent the case when plotted against La, which supports a 
genuine three-body force of Y-type. In Figl^the fits by the y-Ansatz and by the A-Ansatz 
for the quarks in the equilateral triangle are shown in the top and bottom parts, respectively. 
Note that for the quarks in the equilateral triangle these two fits are essentially the same 
with = :^cry^ and equal selfenergy and Coulomb coefficient. 

In Fig0 and FiglHl we show further comparison of our data with A and y-Ansatze for 
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-| I Li I I I I 

0.5 1 1.5 

FIG. 5: The abelian three-quark potential and A and y-Ansatze in fuh QCD for quarks in equi- 
lateral triangles as function of quark separation d = Ly/V^ = The solid line is a fit to the 
data, the dotted line is the A-Ansatz prediction eq. (|2^-i|) and the dashed line is the y-Ansatz 
prediction ea. (|^3|) . For both Ansatze ea. H24() was used. 



full and quenched QCD. The data for the three-quark potential are plotted as a function 
of distance d for equilateral triangle. In the same figures the curves showing respective Y- 
Ansatz and A-Ansatz predictions are plotted. We fix all parameters in the Ansatze using 
relations ()19|) . ()22|) and (j23). We see that for both full and quenched QCD at distances 
(i < 0.5 fm the three quark potential data agree with the A-Ansatz, while at larger distances 
it agrees with F-Ansatz up to an additive constant indicating that the string tension o'y'^^ 
is equal to a'^^ as was already discussed above. Similar findings were presented for the 
quenched nonabelian potential in Thus we conclude that our data for the abelian 
potential confirm the F-Ansatz for large distances. The agreement with the A-Ansatz 
at short distances, which was also observed in P|, is probably a coincidence since the A- 
Ansatz prediction eq. (j22j) is formulated for large quark separations. On the other hand, the 
proximity of the potential to the A-Ansatz at distances which are relevant for the spectrum 
calculations might be important for the phenomenologists since the calculations with the A- 
Ansatz potential are much simpler. The disagreement with the Y-Anzatz at small distances 
was first clearly observed in Ref. P|. One can guess that the finite size of the junction play 
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0.5 1 1.5 



FIG. 6: Same as FigOfor quenched QCD. 



though our data for the static potential at 
5| we are not in a position to make strong 



the role in appearance of this discrepancy. A 
small distances behave similar to that of Ref. 
statements about the short distances since we are using the abelian projection which, as 
many earlier observations suggest, gives correct description of the static potentials at large 
distances only. 

Although the results for the static potential are in favor of the F-Ansatz, the difference 
from the A-Ansatz prediction is rather small. Thus it is worthwhile to study the color 
flux distribution. In Fig|7|we show the distribution of the color electric field E^^ , and its 
surrounding monopole currents k^^ , on the 24^48 lattice in full QCD. The time direction 
of the Wilson loop has been taken in one of the spatial directions of the lattice. Points on 
the hyperplane orthogonal to the time direction of the Wilson loop are marked by {x,y, z). 
The static quarks are placed at (x, y, z) = (20, 10, 8), (25, 18, 8) and (30, 10, 8), respectively, 
i.e. they lie in the {x, y) plane. The color index of the electric field operator (cf. eq. (fTBjl ) is 
identified with the color index of the quark in the bottom-right corner (in the center bottom 
figure). Note that the sum of the electric field over the three color indices vanishes at any 
point. As expected, the flux emanates from the quark in the bottom-right corner and at 
about the center of the 3Q system splits into two parts. The flux lines are schematically 
drawn in FiglHl A similar picture holds for the top and bottom-left quark and their respective 
fluxes. In the adjacent figures we show the monopole current in the planes perpendicular 
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FIG. 7: Distribution of the color electric field E^^ in the (x, y) plane on the 24^ 48 lattice (center 
bottom figure), together with the monopole currents k^*^ in the {x,z) and {y,z) planes (adjacent 
figures), respectively, at the position marked by the respective solid lines. The magnitude of E^'^ 
and k^^ is indicated by the length of the arrows. 



FIG. 8: Schematic view of the color electric field. 



12 



to the electric flux hnes, i.e. the {x, z) and {y, z) planes. They form a solenoidal current, 
as in the case of the QQ system, in agreement with the dual superconductor picture of 
confinement. 

We may decompose the abelian gauge field into a monopole and photon part according 



to the definition |20|,|21| 

^.(s,/i)=r-(s,^) + ^f(s,/i), 

( ^ (26) 
^r"(^, = E.' D{s - s')VVm,{s\ a, ^^) , 

where D{s) = A~^{s) is the lattice Coulomb propagator, is the lattice backward deriva- 
tive, and mi{s, /i, z/) counts the number of Dirac strings piercing the plaquette Ui{s, /i, z/). If 
one computes ki{*s,fi) from 6Y'°'^{s,fi) one recovers almost all monopole currents. In FigE] 
we see that the monopole part is largely responsible for the linear behavior of the potential, 
as was found already in case of the QQ potential The ratio of monopole to abelian 
string tension turns out to be 0.81(3). 

In Figini we show the distribution of the abelian color electric field photon parts. The 
photon part shows a Coulomb-like distribution, while the monopole part has no sources. 
Outside the fiux tube the monopole and photon parts of the color electric field largely cancel. 
The middle figure shows clearly that the fiux lines are attracted to a Y-type geometry. 

In FigEHwe show the action density of the 3Q system in full QCD. Also shown is the 
monopole and photon part of separately. Let us first look at the (full) abelian density. 
It clearly displays a Y-type geometry of the color forces. This is, of course, indistinguishable 
from a geometry of purely two-body forces with strongly attracting flux lines. The monopole 
part of the action density shows no sources. Apart from that, it appears that the action 
density originates almost entirely from the monopole part. The sources show up in the 
photon part of the action density as expected. 

We have done similar calculations (to the ones shown in FigsIH IHl lEI) in quenched 
QCD as well. Part of our flndings have been reported in 0, and we refrain from repeating 
them here. Qualitatively, we found the same results as in full QCD at niT^/mp ^ 0.7. In 
Fig. ^2we compare the action density of full and quenched QCD. We see that at the center 
of the flux tube the action density in full QCD is slightly higher than for the quenched case, 
while the shapes are rather similar. The same feature has been observed for the QQ flux 



tube We have estimated the width 6 of the flux tube using a Gaussian flt 



Ifl. The 
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FIG. 9: Distribution of the abelian color electric field E^^ (top) broken into monopole (middle) 
and photon parts (bottom) on the 24^ 48 lattice in full QCD. The color index of the electric field 
operators corresponds to that of the quark in the bottom right corner. Only part of the lattice is 
shown here. 

result is 5 = 0.30(4) fm and 0.36(11) fm in full and quenched QCD, respectively. This is 
to be compared with the width of the QQ flux tube, which turned out to be 0.29(1) fm 
in full and quenched QCD [10]. We found that the width increases closer to the junction. 
So the numbers quoted above are only to tell that the width of the baryon flux tube, away 
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FIG. 10: The abelian action density of the 3Q system in full QCD, together with the monopole 
and photon part. 
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from the junction is not very different from that of the QQ flux tube. For a more precise 
determination of the width larger quark separation is necessary. 
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FIG. 11: The action density p^^rQ of Fig. plotted across the flux tube at the distance 3a from 
the quark and 2a from the junction. 



It is interesting to compare the action density shown in Fig. with the action density 
constructed out of three QQ flux tube action densities multiplied, in agreement with (|^. 
by a factor ^ to take into account that we are dealing with pairs of quarks rather than with 
QQ pairs. Such a comparison has been done in Ref. ^ for the Potts model. For the QQ 
action density we used the results of Ref. Q] . The resulting density is shown in Fig. [T21 
Figures ITUl and IT^ look rather different. The most important difference is that the measured 
density has a bump in the center, while the A-Ansatz density has a dip. This comparison 
gives further support to the F-Ansatz. 



IV. BARYONIC FLUX AT FINITE TEMPERATURE 

We expect the flux tube to disappear and the color electric field to become Coulomb- or 
Yukawa-like above the finite temperature phase transition and when the string breaks in 
full QCD. This phenomenon has been observed in case of the QQ system in the pure SU(2) 
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FIG. 12: The abelian action density of the 3Q system as predicted by the A-Ansatz in fuh QCD. 



gauge theory for temperatures T > j2j| and in full QCD for T just below and above 
Tc j2j|. Throughout this Section we shall use the Polyakov loop (fT^ to create a baryon. 

In FigElwe show the baryon potential on the 16^8 lattice at /? = 5.2 for several values 
of K. At this (3 value 

T oc exp(-2.81/K) . (27) 

Increasing k thus increases the temperature. We cross the finite temperature phase transition 
at K = 0.1344 We see that the potential flattens off while we approach the transition 
point. However, the distances we were able to probe are not large enough to make any 
statement about string breaking. 

To compute the action density p^^ and the electric field and monopole correlators E^^ 
and fc^'^, respectively, we need to reduce the statistical noise. We do that by averaging over 
time slices and using extended operators 



(28) 



+ pf{s -!-?>)+ pfis -2-3)+pfis~ i) 
+ pf{s-2)+pf{s-3)}, 



E^Q{s, i) ^ ^{E^'^is, i) + E^«(s - m)} , (29) 
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FIG. 13: The monopole part of the baryon potential at finite temperature in full QCD as a function 
of Ly [T < Tc) and La (T > Tc), respectively. 



eQ{*s,^i) -{k'^^Cs^fi) + - 2),^)} , /i = 1,3 (30) 

where (again) we have assumed that the quarks he in the {x, y) plane, and we consider the 
monopole current in the (z, x) plane. 

In Fig^^we plot the abelian action density in the deconfined phase at k = 0.1360. As 
was to be expected, the action density shows three Coulomb-like peaks at the position of 
the quarks, similar to the photon part of the action density at zero temperature as shown 
in FigHH 

In Fig El we show the monopole part of the electric field, averaged over the color compo- 
nents, and the accompanying monopole current for three values of k, corresponding (from 
left to right) to the confined case, to T ^ Tc and to the deconfined phase. In the confine- 
ment phase (k = 0.1335) we find the fiux to be of Y-shape, similar to the zero temperature 
case where we used Wilson loop correlators. Note that the Polyakov lines do not have a 
Y-shape junction like the Wilson loop does, which excludes the possibility that the fiux is 
being induced by the color lines. Just below [k = 0.1344) we still see a Y-shape fiux, 
while in the deconfined phase {k = 0.1360) the electric field becomes Coulomb-like. 
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FIG. 14: The Abelian action density in the deconfined phase at k = 0.1360. 
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FIG. 15: The color electric field (top) and monopole currents (bottom) on the 16^ 8 lattice at 
K = 0.1335 (left), 0.1344 (middle) and 0.1360 (right). The three quarks lie in (what we call) the 
(x,y) plane. The bottom figures show the monopole currents in the (x, z) plane at the position 
marked by the solid lines in the top figures. 
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V. CONCLUSIONS 



We have studied the 3Q system in the maximally abelian gauge in full QCD at zero and at 
finite temperature. Among the quantities we have looked at are the abelian baryon potential 
as well as the fiux distribution and the action density. While on the basis of the potential 
it is hard to decide whether the long-range potential is of A- or Y-type, the distribution of 
the color electric field and the action density clearly shows a Y-shape geometry. As in the 
QQ system, we identified the solenoidal monopole current to be responsible for squeezing 
the color electric fiux into a narrow tube. Little difference to the quenched theory was 
found. In the deconfined phase the fiux tube disappears, and the color electric field assumes 
a Coulomb-like form. Our results are in qualitative agreement with the predictions of the 
dual Ginzburg-Landau model 



24| : the baryon fiux has F-shape, and the solenoidal monopole 



currents are clearly observed. 
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